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(57) A video decompression system includes block 
error concealment means (22,41,42) for generating 
substitute data for lost blocks of pixel data. An interpola- 
tor generates both spatially (44) and temporally (40) 
interpolated or predicted blocks of data for concealing a 
lost block. Pixel data in the respective blocks of spatially 
and interpolated data are transformed (43) to coeffi- 
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cients representing frequency spectra. A substitute 
block of coefficients is assembled from coefficients from 
both the transformed blocks, according to a predeter- 
mined criteria The substitute block is transformed back 
to the spatial domain for substitution of lost pixel values. 
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Description 

This invention relates to error concealment apparatus for replacing damaged or lost two dimensional blocks of pixel 
values in reproduced images. 

5 In video signal transmission systems of the type which compress video data on a block basis as for example sys- 

tems which utilize discrete cosine transforms, data loss and/or data damage frequently occurs on a block basis (e.g. a 
16X16 array of pixels). These transmission systems normally do not support enough overhead to correct allerrors, but 
rather depend on error concealment techniques to provide substitute data which closely approximates the lost data. 
Werner Keesen et. al. in United States Patent No. 4,807,033 issued 21 February 1989, disclose an error concealment 

10 system for replacing two dimensional blocks of image data. This system includes apparatus for detecting blocks of lost 
data and replacing the lost image data in the spatial domain with pixel values generated from surrounding pixels by 
interpolation. The surrounding pixels are examined to determine the existence of image edges (gradients) and thereaf- 
ter interpolation is performed to generate a block of pixel values to replace the lost or damaged block of pixel values. 
Directional interpolation, as used herein, is the process of generating pixel values by performing interpolation which 

is favors a particular direction. Refer to FIGURE 1 where the respective squares or blocks represent respective image 
areas forming the basis of block based compression. The dots within the respective squares or blocks represent pixel 
values, in this instance representing white or black image points. The center block which is void of dots represents a 
damaged block for which substitute pixels are to be generated. If for example, pairs of pixels within the two rows of pix- 
els circumscribing the central block are appropriately examined, existing image gradients may be determined as for 

20 example the diagonal white line passing through blocks B1 and B9. Then, based on the gradient information, interpola- 
tion is performed in the direction of the dominant detected gradients. In the Figure, the pixels which fall on a given arrow 
are generated from pixels in the adjacent blocks and occurring along extensions of the respective arrows. 

The Keesen et. al apparatus detects image gradient information and for single gradients performs one directional 
interpolation. If more than one strong gradient is detected, a directional interpolation is performed for each of the strong 

25 gradients and the arithmetic sum of the respective directionally interpolated blocks of pixels is utilized as the substitute 
block. 

Summing the respective interpolated blocks of pixels gives respective pixel values equal weight and tends to blur 
image features. Consider a gray image containing a horizontal white line behind a vertical black line. Vertical interpola- 
tion will generate a block of gray pixels including the vertical black line. Horizontal interpolation will generate a block of 
30 gray pixels including the horizontal white line The arithmetic average of these interpolated blocks will be a block of gray 
pixels including a lighter than black vertical line and a darker than white horizontal line, with a gray area at the intersec- 
tion of the lines. 

The present invention is an improvement over the Keesen et. al. error concealment apparatus and includes block 
error concealment means for generating substitute data for lost blocks. An interpolation apparatus generates both spa- 

35 tially and temporally interpolated or predicted blocks of data for concealing a lost block. Pixel data in the respective 
blocks of spatially and interpolated data are transformed to coefficients representing frequency spectra. A substitute 
block of coefficients is assembled from coefficients from both the transformed blocks, according to a predetermined cri- 
teria The substitute block is transformed back to the spatial domain for substitution of lost pixel values. 

In a particular embodiment, the relative correlation of parameters in blocks adjacent a missing block is determined. 

40 if the relative correlation falls within a predetermined range, the foregoing substitute block is utilized for block error con- 
cealment. On the other hand, if the relative correlation exceeds or is less than the range, then the temporally predicted 
block and the spatially interpolated block are utilized for block error concealment respectively. 

BRIEF DESCRIPTION OF THE DRAWINGS 

45 

FIGURE 1 is a pictorial representation of a portion of the blocks of pixel values representing an image, which rep- 
resentation is useful in describing the invention. 

FIGURE 2 is a block diagram of apparatus for decompressing video signal and including block based error conceal- 
ment apparatus embodying the present invention. 
so FIGURE 3 is a pictorial representation of the process of error concealment use in the present invention 

FIGURE 4 is block diagram of exemplary concealment apparatus which may be implemented for element 26 of 
FIGURE 2. 

FIGURE 5 is a pictorial representation of transform coefficients useful in describing a method for determining dom- 
inant image edges or direction, and describing directional spatial interpolation. 
55 FIGURE 6 is a pictorial diagram indicating the location of blocks from which motion vectors are accessed for eval- 
uation of image motion. 

FIGURE 7 is a flow chart showing the operation of the apparatus of FIGURE 4. 
FIGURE 8 is a flow chart, illustrating in part, the process 61 of FIGURE 7. 



2 



BNSDOCID: <EP 072791 0A2J__> 



EP0 727 910 A2 



FIGURE 9 is a block diagram of apparatus for generating combined spatio-temporal interpolated error concealment 
blocks. 

FIGURE 10 is a block diagram of apparatus for calculating a variable parameter "e" used in the FIGURE 9 appara- 
tus. 

In order to perform block based error concealment it is necessary to be able to determine lost or damaged blocks 
of signal. To this end it is assumed that the transmitted signal to be processed in the exemplary embodiments is format- 
ted in transport packets. Each packet includes a header, a payload, and continuity and/or check bytes. The header 
includes sufficient information to identify the spatial locations, in the reproduced image, of the respective blocks of com- 
pressed video data carried in the payload. The check bytes may be of a form to permit error detection of substantially 
io all potential errors and at least partial error correction. An example of this type of signal format may be found in United 
states Patent No. 5,247,363. Of course there are other known methods for detecting lost blocks which may also be 
implemented. 

Refer to FIGURE 2 which illustrates a general form of a receiver for processing the above described signals. Trans- 
mitted signals, from for example, an antenna (not shown) are applied to a tuner/demodulator 10, which provides base- 
15 band signal. The baseband signal is applied to a forward error correction, FEC, decoder 12, which examines the 
transmitted signal and corrects errors incurred in the transmission channel, according to its error correction capability. 
Error corrected data from the FEC decoder 1 2 is applied to a rate buffer 1 4 which accepts data at a relatively constant 
rate from the FEC decoder 12 and outputs data on demand to the subsequent processing elements. Data from the 
buffer is coupled to a frame check sequence, FCS, decoder 16, which examines the error corrected data for uncorrected 
20 errors according to FCS check bits (or continuity check count values, for example) appended to the transport packets. 

The FCS decoder 16 passes transport packets, D, to element 18 along with an error signal, ED, indicating whether 
respective transport packets contain data errors. Element 18, responsive to the error data, discards transport packets 
which include data errors. Error free transport packets are separated into video data and transport header data. The 
transport header data and error data, ED, are utilized to determine what data is lost. For more specific information 
25 regarding these processes the reader is invited to examine U.S. Patent application (RCA 86409A) Sn. 08/01 7,455, filed 
•February 12, 1993. and entitled "APPARATUS FOR CONCEALING ERRORS IN A DIGITAL VIDEO PROCESSING 
SYSTEM". Element 18 passes the error free video data and the error tokens ED to the decompressor 20. The error 
tokens may provide several functions. For example, the particular error token may condition the decompressor to simply 
skip over lost blocks in particular frame types. In this instance, decompressed data from the prior frame is simply 
30 repeated from the collocated image area, i.e., a form of temporal interpolation or block replacement. Alternatively, other 
error tokens are passed to a further error concealment processor 26. 

The decompressor 20 utilizes frames of memory 22 in the decompression process. This memory may be arranged 
to store frames of both compressed and decoded video signal for use in the process of motion compensated predictive 
decoding, and error concealment Decoded image information is subsequently transferred to a display memory 24 from 
35 which it may be accessed in raster format for display. The decompressor in this arrangement is conditioned to store 
decoded motion vectors, for at least a portion of a frame, in the memory 22 for use in error concealment. 

The error concealment element 26 interacts with the data in memory 22, and generates substitute image data for 
blocks of pixel values identified by the error tokens/Luminance and chrominance signal components are decompressed 
separately but by similar processes and then recombined for display. The following discussion regarding error conceal- 
40 ment is directed to concealment of errors in the luminance component. Errors in the chrominance components may be 
concealed in similar fashion by similar parallel apparatus not shown. However, parameter evaluation may be performed 
on a single signal component (e.g. luminance) and the results of the evaluation used for determining the type of block 
error concealment for all components. 

FIGURE 3 shows pictorially the functions of one embodiment of the error concealment apparatus 26 of FIGURE 2. 
45 In FIGURE 3, blocks 30 and 31 represent memory from which decompressed data may be accessed. The data from 
element 30 is data which surrounds the lost block B. The data, TR, from element 31 , preferably is a temporal prediction 
of the lost block using motion vectors from the block vertically above the lost block to determine the predicted values. 
Alternatively, block TR may represent data from a prior frame collocated with block B. It should be realized that the 
blocks to the right and left of block B may also be lost or missing if they were contained in the same transport block. 
so However, even if this were the case, the block to the left of block B may be a replacement block if the system performs 
error concealment from left to right, in which case only the block to the right of block B will also be missing. The data 
from elements 30 and 31 are evaluated to determine the mode of error concealment to be executed for the missing 
block B. In order to facilitate this evaluation, it may be beneficial to first replace the missing data of block B and possibly 
the blocks to the right and left of block B (if they too are missing) with some predetermined value, representing gray for 
55 example, or possibly with the average brightness level of the block TR. 

Two or three types of image evaluation are executed 32, to determine the concealment mode to be performed. A 
first evaluation is for image motion. An example of image motion evaluation, to be described below, looks at the magni- 
tude and relative correlation of the motion vectors for the blocks above and below the missing block. A second evalua- 
tion is of dominant image features such as edges or gradients. An example of edge evaluation, to be described below, 
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D£ rermines dominant edges in the area surrounding the lost block B, and dominant edges in the collocated block TR, 
ena also examines the correlation between the edges in blocks B and TR. A third evaluation takes a general measure 
c race differences between an area of the current frame being decoded and a collocated area of a prior frame. From 
trv foregoing evaluations a closeness measure is derived, and depending upon the relative value of this measure one 

s of three modes of error concealment is performed. 

The first mode of concealment, 35, executed when high correlation exists between the current and the prior frame 
is simple temporal replacement of the lost block B with a temporally predicted block, TR. A second mode of replace- 
ment, 34, executed when low correlation is indicated by the closeness measure, is concealment by spatial interpolation. 
The interpolation in this instance is directional, according to the direction of the dominant edge or edges determined in 

10 the neighboring image areas. A third mode of replacement, 33 combines the prior two modes. The combining is per- 
formed in the spatial frequency domain wherein frequency components of the dominant edge are non-additively com- 
bined with frequency components of the block TR exclusive of the frequency components representing the edge. These 
combined frequency components are inverse transformed to generate concealment data for the lost block. 

This type of error concealment process permits of a very advantageous trade off between generating high resolu- 

75 tion replacement blocks and blocks with reduced artifacts in the presence of motion, through the range of images with 
no motion and images with significant motion. The combination of spatial and temporal interpolation for images having 
intermediate correlation, advantageously provides foreground images with lesser resolution but also with lesser motion 
artifacts and background images with high spatial resolution tending to give the entire substitute image block the 
appearance of high spatial resolution. 

20 FIGURE 4 illustrates, in block diagram form, exemplary error concealment apparatus. The illustrated functions may 

-*e programmed into the software of a special purpose processor such as the CLM41 00 Multimedia Accelerator availa- 
oie from C-Cube Microsystems, Milpitas, CA 95035, or they may be realized by individual circuit elements as indicated. 
\ ne apparatus includes element 43 designated as a transform, inverse transform device. The transform, and inverse 
transform will be described in terms of a Fast Fourier Transform, FFT, and Inverse Fast Fourier Transform, IFFT respec- 

25 tively. However it will be appreciated by those skilled in the art of signal processing that other transforms may be used, 
in particular the Discrete Cosine Transform, DCT and its inverse, IDCT for example. 

In FIGURE 4 error data from the decompressor 20 triggers a concealment controller 40, which may be a microproc- 
essor, to generate interpolated data for a lost block of data. In furtherance of this process, controller 40 accesses motion 
vectors and decompressed image data which surrounds the lost block from the memory 22, and collocated data from 

30 the prior image frame corresponding to block TR. This data may be transferred to a working memory 41 , to preclude, 
at least in part, memory access contention in memory 22 during the concealment processing operations. 

The image analysis is performed by the digital signal processor 42 under control of the concealment controller 40. 
The digital signal processor may be of the type TMS320 series of DSP's available from Texas Instruments Inc., Dallas, 
Texas, for example. Appropriate image data is accessed from the working memory to perform the edge or gradient anal- 

35 ysis for the image area surrounding the lost block. For example, if the lost block contains a matrix of N by N pixels, a 
superblock of 2N by 2N pixels, with the lost block in its center, may be accessed from memory 41 . The lost pixel data 
may be substituted with a mid gray value prior to analysis. The dominant gradient of the superblock is determined, and 
its angle, Theta(M) is recorded. The predicted block (TR) of pixel values from, for example, the prior frame or closest 
anchor frame, is accessed from the memory 41 and applied to the digital signal processor 42 for gradient analysis. The 

40 dominant gradient of the block (TR) is determined, and its angle, Theta(T) is recorded. The values Theta(M) and 
Theta(T) are then applied to the digital signal processor for correlation according to the function; 

r d = cos(Theta(M) - Theta(T)) 

45 wnere r d is designated an image gradient correlation measure. 

In the spatial domain, the DSP 42 is programmed to determine significant image gradients represented by the pix- 
els surrounding the lost block, and then to select the top one or two gradient directions that tend to identify the predom- 
inant local image contours (or most strongly characterize image edge orientations). 

Image edges or gradients may be determined as follows.. The local edge gradients g x and g y for a respective pixel 
so location x(t.j), are computed from the equations: 

9x = x Mj+1 " X MJ-1 + x ij+1 " x iJ-1 + x i+1j+1 " x i+1j-1 ( 1 ) 
9y = X <+1,j-1 " x i-1j-1 + x i+1,j " x i-1.j +x i+1 j+1 " x i-1,j+1 (2) 

This is equivalent to applying the following 3X3 Prewitt convolutional mask operators: 
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The magnitude, G, and angular direction, 4>, of the gradient at (i j) are given by: 

G=*/9x 2 +9y 2 *=tan- 1 (g y /g x ) (3) 

These gradient measures are computed for each (i,j) coordinate in the neighborhood surrounding the missing or dam- 
aged block. The value of the respective gradient angles are rounded to the nearest 22.5 degrees, for example, and thus 
corresponds to one of eight directional categories D1 - D8. If a line drawn through a respective (i.j) neighborhood pixel 
coordinate with the direction determined by the gradient angle will pass through the missing block, that particular (ij) 
pixel coordinate is given a vote. The directional categories D1 - D8 of the voting coordinates are mutually exclusively 
accumulated by summing the magnitudes of the associated directional gradients. That is, the magnitudes of the voting 
coordinates having directional category D1 are accumulated to form one sum, the magnitudes of the voting coordinates 
having directional category D2 are accumulated to form a second sum and so on. The directional categories exhibiting 
the largest accumulated sums determine pixel interpolation directions. 

An alternative method of detecting gradient information will be described with reference to FIGURE 5. The circles 
in FIGURE 5 represent FFT coefficient values for a superblock, arranged in relation to a vertical frequency axis. f v , and 
a horizontal frequency axis, f H . The coefficients represent spatial frequency coefficients for an image and in this exam- 
ple correspond to a 16 by 16 pixel matrix. The coefficient at the intersection of the f v and f H axes corresponds to DC or 
the average brightness level over the matrix. Coefficients to the right of the vertical axis f v represent increasing horizon- 
tal frequencies and coefficients above the horizontal axis f H represent increasing vertical frequencies. Coefficients to 
the left of the vertical axis and below the horizontal axis represent the conjugates of the foregoing frequency compo- 
nents. 

In FIGURE 5 exemplary coefficient groupings designated 135° FILTER. 90° FILTER and H FILTER are shown. The 
coefficients in the group designated H FILTER include selected coefficients which correspond generally to horizontal 
image gradients. The coefficients in the group designated 90° FILTER include selected coefficients which correspond 
generally to vertical image gradients. The coefficients in the group designated 135° FILTER include selected coeffi- 
cients which correspond generally to image gradients at 135 degrees to the horizontal etc. Therefore, it may be seen 
that respective groupings can be arranged for a plurality of image gradient directions. To determine the dominant image 
gradient direction, one need only sum the coefficients in the respective groupings (each of which represents a particular 
direction), normalize the sums and compare the respective normalized sums to determine the largest normalized sum. 
The angle Theta(M) associated with the group having the largest normalized sum is the dominant image gradient. 

Gradient analysis may be performed on the superblock by applying the appropriate pixel data from memory 41 to 
the transform apparatus 43 to generate the corresponding coefficients which will be stored in the memory 41 for the 
above described analysis, which may be performed in the digital signal processor (DSP) 42. The image gradients of the 
predicted block TR may be similarly determined by transforming the pixel data of block TR in element 43, storing the 
coefficients in memory 41 and analyzing respective groupings of coefficients. If this method of gradient analysis is uti- 
lized, the transform coefficients will be retained in memory for further processing described below. 

Image motion analysis is also performed. An exemplary motion analysis generates six motion vector correlation 
measures and forms the weighted average of the six as a directional magnitude correlation measure r m . The motion 
vectors utilized in the analysis are the vectors associated with the blocks designated TOP-1 , TOP, TOP+1 , BOT-1 , BOT. 
and BOT+1, located above and below the missing block as shown in FIGURE 6. Note, for MPEG data, motion vectors 
are associated with macroblocks which include a matrix of 2 by 2 luminance representative blocks each representing a 
matrix of 8 by 8 pixels. If block replacement (concealment) is performed on 8 by 8 pixel units, one directional magnitude 
correlation will serve four blocks. However, since in the MPEG environment it is unlikely that only one block of a mac- 
roblock will be lost, but rather that all blocks of respective macroblocks will be concurrently lost, all analysis will be made 
on a macroblock basis and concealment will be performed on a macroblock basis. It should be assumed that for an 
MPEG environment that the blocks drawn in FIGURE 6 represent macroblocks. 

The motion vectors associated with the blocks neighboring the lost block are accessed from the memory 41 and 
applied to the DSP 42. The DSP 42 is conditioned by the controller 40 to calculate the following correlation factors: 

V TOP* V TOP-1 

fl = IV topIIVtop-iI 
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where V (i) is the motion vector associated with block i. For example V B oi-i «s the motion vector associated with block 
20 BOT-1 . The directional magnitude correlation measure r m is generated according to the relation; 

r m =(K1)Max(r 5 , 0) + (K2)Max(Max(r v r 2 ), 0) + (K3)Max(Max(r v r 2 ), 0) + (K4)r 6 . 

Exemplary values for constants K1 -K4 are 0.6, 0. 1 5, 0. 15 and 0. 1 respectively. Note that to lessen computations r m may 
25 be calculated, with only slight loss in effectiveness, using only r 1 , r 2 , and r 5 with appropriate K scale factors. 

The foregoing vector analysis assumes that respective blocks have motion vectors of similar form. However the 
MPEG compressed video signal protocol supports blocks without motion vectors (I frames), blocks with forward motion 
vectors (P frames) and blocks with forward and backward motion vectors (B frames), and frames with admixtures of 
each type of motion vector. The foregoing motion vector correlation analysis only makes sense if the respective pairs of 
30 macroblocks having motion vectors included in each of the equations r A -r 6 are of the same type. Therefor a test is per- 
formed on the respective pairs of motion vectors before inclusion in the equations. The tests may be performed by the 
controller 40 or the DSP 42 or a combination of both. 

An alternative motion vector correlation calculation, considering the possibilities of different types of motion vectors, 
is illustrated in FIGURE 8. In FIGURE 8, the block types of pairs of the blocks illustrated in FIGURE 6 are accessed 
35 {80}. The blocks are tested {81} to determine if they include the same type of vectors. If they do not, the vector correla- 
tion v(n) for the pair of blocks is set to zero {83}. If the blocks contain the same type of vectors, they are tested {82} to 
determine if they are both I blocks. If they are, they do not contain vectors and the vector correlation v(n) for that pair of 
blocks is set to zero {84}. If the pair of blocks are not both I blocks {82}, they are tested {85} to determine whether they 
both contain backward vectors. If they do, a vector correlation v(n) is calculated {86}according to the equation 



40 



55 



v(n) = cos(angle between vb(i) and vb(i+1)) 



where vb(i) and vb(i+1) correspond to the backward vectors from blocks (i) and (i+1) respectively. 

If at {85} the blocks do not both contain backward vectors, a test is made {87} to determine if both blocks contain 
45 only forward vectors. If they do, a vector correlation v(n) is calculated {88} according to the equation; 

v(n) = cos(angle between vf(i) and vf(i+1)) 

where vf(i) and vf(i+1) correspond to the forward vectors from blocks (i) and (i+1) respectively. 
so If both blocks do not contain only forward or only backward vectors at this juncture they must both contain both for- 

ward and backward vectors. A vector correlation v1 is calculated for the backward pair of vectors and a vector correla- 
tion v2 is calculated for the pair of forward vectors, where 



v1 = cos(angle between vb(i) and vb(i+1)); and 
v2 = cos(angle between vf(i) and vf(i+1)). 

The vector correlation v(n) for the pair of blocks is the larger of v1 and v2 {89}. After each pair of blocks ( for exam- 
ple TOP, TOP-1 ; TOP, TOP+1 ; BOT, BOT-1 ; BOX BOT+1 ; BOT, TOP as illustrated in FIGURE 6) is examined and a cor- 
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responding vector correlation is determined, a check is made {90} to determine if all desired pairs of blocks (or 
macroblocks) have been processed. If they have not, the next pair is accessed {80} and processed. If they have all been 
processed a correlation measure r m is generated according to the equation ; 



N 

r m = £K(n)v(n) 



where K(n) for forward and backward types are the same value, but K(n) for a v(n) associated with a pair of blocks con- 
io taining both forward and backward vectors is larger. 

The temporal correlation measure or error analysis, r p , is performed according to the equation; 

15 15 15 15 
( X I TOPijTOPPij) JroPij ^TOPPij 
j=o i=0 i=0 j=0 



15 



To = 



STD ( TOP ) STD ( TOPP ) 



20 



where TOPjj refers to the pixel values in the block above the missing block assuming a block of 16 by 16 pixels, TOPPjj 
refers to the pixel values for a predicted block corresponding to the block above the missing block, and STD stands for 
standard deviation. In determining the predicted block TOPP. the motion vectors from block TOP are used unless the 

25 block TOP is an I type block, in which case a zero valued motion vector is used. The calculation for r p may be performed 
in either the DSP 42, the controller 40. or a combination of both. 

The same calculation may be made for the lower block BOT and its temporal prediction BOTP or only the first cal- 
culation may be performed as desired. If a calculation is performed for both the top (TOP) and bottom (BOT) blocks the 
larger of the r p values is selected as the temporal correlation measure r p . 

30 The different correlation measures r p , r m , and r d are added to generate a general correlation parameter r, according 
to the relation; 

r-PKr p ) + p2(r d ) + p3(r m ) 

35 where p1 , p2, and p3 are constants having exemplary values of 0.6, 0.2 and 0.2 respectively. Error concealment mode 
is determined by the relative magnitude of r. If r is greater than a threshold TH1 , indicating a high degree of neighboring 
image correlation, concealment of the missing block is by temporal replacement using a predicted block from a prior 
frame. If r is less than a threshold TH2, (TH2 < TH1) indicating relatively low neighboring image correlation, conceal- 
ment of the missing block is by substituting a block generated by spatial interpolation. If TH2 ^ r ^ TH1 concealment is 

40 by a combination of spatial interpolation and temporal replacement to be discussed below. Exemplary values of TH1 
and TH2 are 0.9 and 0.4 respectively. 

It is not necessary to calculate and utilize all three of the correlation measures and acceptable results may still be 
achieved using for example only two of the correlation measures. Alternatively other correlation values/measures may 
be added to the mix forming the value r. In addition one or the other of the coefficients pi , p2 and p3 may be set to zero. 

45 The foregoing processes are reviewed in the flowchart of FIGURE 7. It will be appreciated that the order of deter- 
mining the coefficient measures r m , r p , and r d is interchangeable. In FIGURE 7 the measure r p is first calculated {60}. 
Then the appropriate motion vectors are accessed and the correlation measure r m is determined {61}. The pixel values 
for blocks neighboring a lost block are accessed {62} and the edge gradients r d are determined {63}. The respective 
measures r p , r m and r d are weighted and summed to generate the correlation measure of the spatially and temporally 

so surrounding blocks {64}. 

The value of the correlation parameter r is compared to two threshold values TH1 and TH2 in the concealment con- 
troller 40. If r is greater than TH1 a substitute block for the bad or missing block is generated by temporal replacement. 
This is performed in the controller 40, by substituting the block TR for the missing block in the display portion of memory 
22, for the current frame period. If the value r is less than TH2. then concealment is by substituting a spatially interpo- 

55 lated block. This may be performed by accessing the superblock from the working memory 41 , and applying this pixel 
matrix to the spatial interpolator 44. The interpolator 44 may generate a substitute block by known interpolation meth- 
ods, including directional or bi-directional spatial interpolation in accordance with the directions of the dominant image 
gradient or gradients. Alternatively, if TH2 =s r ^ TH1 , concealment is accomplished with a combination of spatial and 
temporal interpolation, which will be described with reference to FIGURE 9. 
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In FIGURE 9 the block 900 corresponds to the superblock previously discussed for generating the dominant edges 
or gradients. Element 901 represents apparatus for detecting the dominant edges or gradients. Indicia indicating the 
direction of the dominant edges are applied to a gradient mask 904. The superblock is applied to an FFT 903, and the 
transform coefficients generated by the FFT are applied to the gradient mask 904. Mask 904, responsive to the indicia 

5 from the detector 901 , selects only those coefficients associated with the dominant edges. Referring to FIGURE 5, if 
the dominant edge is at 135 degrees, the gradient mask 904 will form a superblock coefficient matrix comprising those 
coefficients which fall within the group designated 135° FILTER (in FIGURE 5), and with the remaining coefficients set 
to zero values. This matrix is applied to an inverse FFT 906 and the coefficients are transformed to generate a super- 
block spatial pixel matrix 908. (Item 908 represents a portion of working memory 41 , FFT 903 and IFFT 906 correspond 

to to element 43 in FIGURE 4 and the gradient mask may be realized by appropriate programming of controller 40.) The 
superblock is applied to a selector 910 which extracts the matrix of pixels corresponding to the missing block. That is, 
the pixel values of the superblock which form part of the neighboring blocks are eliminated leaving only the missing 
block M. This block M corresponds to a directional ly filtered block and may be utilized as the spatially interpolated sub- 
stitute block when such substitution mode is signaled by the parameter r. For the current purposes, the block M is 

75 applied to the FFT 914 which generates transformed block M* (916). In addition the predicted block TR from a prior 
frame is accessed from memory 912, and applied to the FFT 914 which generates the transformed block TR' (918). 
Respective FFT coefficients of blocks M' and TR' are applied to a subtracter 920 which forms a matrix of the absolute 
values of the differences of the respective corresponding coefficients. The respective differences |M'[i,j] - TR'[i,j]| are 
compared against a variable "e M in a comparator 922. The output of the comparator 922 controls a selector 926 which 

20 forms a matrix of coefficients Y[ij] (927) from the coefficients in block M' and TR'. If the difference corresponding to 
coefficient [i,j] is less than or equal to "e", the selector 926 is conditioned to select the M[i,j] coefficient in block M' for 
use as the Y[i,j] coefficient. Alternatively, if the difference corresponding to coefficient [i.j] is greater than V, the selector 
926 is conditioned to select the TR[i,j] coefficient in block TR' for use as the respective Y[i.j] coefficient. The Y[i,j] matrix 
is applied to an inverse FFT element 928 which generates the substitute block of pixel values. This matrix may be 

25 returned to the memory 900 and the process iterated if desired. (Note, the elements 910 and 926 may be subsumed in 
element 40 of FIGURE 4, elements 920 and 922 may be subsumed in element 42, and elements 914 and 928 may be 
subsumed in element 43.) 

The dashed box in the block 927 is meant to indicate that the coefficients of the block Y[i,j] may be low pass filtered. 
Low pass filtering may be accomplished simply by ignoring or excising the higher frequency coefficients in the respec- 
30 tive blocks. The low pass filtering may be performed as part of the selection process in the selector 926 and may be 
made dependent upon the value of the correlation measure r. Low pass filtering as a function of r advantageously less- 
ens potential artifacts arising from images having lesser spatial and temporal correlation. Larger values of r would 
impose little or no filtering and smaller values of r would impose greater filtering. Alternatively the filtering may be per- 
formed at the FFT 914, to lessen the amount of processing required of elements 920 and 922. In this instance the 
35 blocks M* and TR' will be independently but similarly filtered. 

The selector 926 may also be arranged to perform the selection of the respective blocks of data dependent upon 
the value of r. That is, the selector may include a comparator to compare r against the threshold values TH1 and TH2, 
and select blocks M* or TR* or a combination of M' and TR' according to the relative magnitude of r. 

The variable "e" may be a preprogrammed constant or a value selected by the user. Alternatively, it may be gener- 
ic ated automatically as a function of the signal. One mode of generating the variable "e" automatically is illustrated in FIG- 
URE 10. In FIGURE 10 the respective pixel values of the current block TOP (1001) and the predicted top block TOPP 
(1002) are applied to a subtracter 1003 which generates a matrix of difference values. These difference values are 
transformed in a two dimensional FFT 1 004 and applied to element 1005 which generates the variable "e" as a function 
of the differences ( here as a function of the transform coefficients of the difference values F(DIFF). An exemplary func- 
45 tion F(DIFF), is to set variable "e" equal to the 90% value of the transform coefficients of the difference values. Alterna- 
tively variable e may be set equal to the value of one standard deviation a of the coefficients weighted by a scale factor 
6 which may be determined empirically via experiment. 

e = 6 times a 

50 

Another value for F(DIFF) = e may be the weighted value of the mean, m, of the coefficients. 

e = K4 times m 

55 where K4 is a constant determined by experiment. Using values of "e" which are a function of the video signal for deter- 
mining the selection of the respective Y[i j] coefficients in the substitute block has the advantage of minimizing errors 
where there are large correlative differences between correlative changes in the image 

An alternative mode of generating the temporal-spatial interpolated substitute block is suggested in part by the 
dashed line from the gradient mask to the selector 926. In this mode the selector 926 is conditioned by the dominant 
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gradient to form a matrix Y[i,j] by selecting from block M\ only the M[i,j] coefficients corresponding to those representing 
the dominant gradient, and selecting coefficients from block TR* to provide the remaining coefficient values. 

Claims 

5 

1 . Apparatus for generating substitute blocks of pixel data for missing or corrupted blocks of pixel data, characterized 
by: 

a source (20) of pixel data representing respective images and a source of error data indicating blocks of 
said pixel data that may be missing or corrupted (hereinafter designated bad pixel data); 
70 a controller (40), responsive to said error data, for accessing from said source, pixel data located spatially 

adjacent to a block of bad pixel data, and a temporally predicted block of data corresponding to said block of bad 
pixel data; 

an interpolator (44) operative upon said pixel data located spatially adjacent to a block of bad pixel data, for 
generating a spatially interpolated block of pixel data; 
75 transform circuitry (43) for providing a matrix of transform coefficients representing said spatially interpo- 

lated block of pixel data, and for generating a further matrix of transform coefficients representing said temporally 
predicted block of data; 

a selector (926), responsive to a predetermined criteria, for forming a still further matrix of transform coeffi- 
cients, by selecting respective coefficients in said still further matrix exclusively from said matrix of transform coef- 
20 ficients representing said spatially interpolated block, or from said further matrix of transform coefficients in 
accordance with said criteria; and 

inverse transform circuitry (928) for inverse transforming said still further matrix to generate said substitute 
block of pixel data. 

25 2. The apparatus set forth in claim 1 characterized in that said interpolator includes digital signal processing appara- 
tus operative upon said pixel data located spatially adjacent to a block of bad pixel data for detecting a direction of 
a dominant image gradient within an image area represented by pixel data located spatially adjacent to a block of 
bad pixel data. 

30 3. The apparatus set forth in claim 2 characterized in that said said selector for forming a still further matrix, generates 
a matrix of transform coefficients including transform coefficients which represent directional image information 
corresponding to said dominant image gradient from said matrix of transform coefficients representing said spa- 
tially interpolated block, and selects transform coefficients from said further matrix of coefficients exclusive of coef- 
ficients which represent directional information corresponding to said dominant edge gradient. 

35 

4. The apparatus set forth in claim 1 characterized in that said said selector comprises; 

circuitry for determining respective differences between corresponding transform coefficients in said matrix 
of transform coefficients representing said spatially interpolated block of pixel data, and in said further matrix of 
transform coefficients; 

40 a comparator for comparing respective differences to a predetermined variable "e"; and 

wherein said selector selects a corresponding coefficient from said matrix of transform coefficients repre- 
senting said spatially interpolated block of pixel data, when the corresponding difference is less than said variable 
"e'\ and selects a corresponding coefficient from said further matrix of transform coefficients when the correspond- 
ing difference is greater than the variable "e". 

45 

5. The apparatus set forth in claim 4 characterized in that said said selector further includes means for generating said 
variable "e" from image information located adjacent said bad data. 

6. The apparatus set forth in claim 4 further characterized by: 

so apparatus responsive to data representing image information of image areas adjacent said block of bad 

pixel data, for determining a measure, r, of correlation of image features in said image areas adjacent said block of 
bad pixel data; 

a source of threshold value TH1 ; 

a comparator for comparing said correlation measure r against the threshold value TH1 and conditioning 
55 said selector to select coefficients only from said further matrix rf the correlation measure, r, is greater than TH1 . 
and conditioning said selector to select coefficients according to said predetermined criteria if the correlation meas- 
ure, r is less than TH1. 
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7. The apparatus set forth in claim 6 further characterized by: 

a source of a further threshold value TH2 where TH2 is less than TH1 ; and 

wherein said comparator is arranged to compare said correlation measure, r, against said further threshold 
value TH2, and to condition said selector to select coefficients only from said matrix of transform coefficients rep- 
resenting said spatially interpolated block of pixel data if the correlation measure, r, is less than TH2, and condition- 
ing said selector to select coefficients according to said predetermined criteria if TH2 < r < TH1 

8. The apparatus set forth in claim 1 further characterized by; 

apparatus responsive to data representing image information of image areas adjacent said block of bad 
pixel data, for determining a measure, r, of correlation of image features in said image areas adjacent said block of 
bad pixel data; 

a source of threshold value TH1 ; 

means for producing a block of substitute pixel data; 

a comparator for comparing said correlation measure r against the threshold value TH1 and conditioning 
said means for producing a block of substitute pixel data to provide said temporally predicted block of data, if the 
correlation measure, r, is greater than TH1 , and to provide a block of data from said selector if the correlation meas- 
ure, r is less than TH1. 

9. The apparatus set forth in claim 8 further characterized by; 

a source of a threshold value TH2 where TH2 is less than TH1; and 

wherein said comparator is arranged for comparing said correlation measure, r, against said further thresh- 
old value TH2, and conditioning said means for producing a block of substitute pixel data to provide said spatially 
interpolated block of pixel data if the correlation measure, r, is less than TH2, and to provide a block of data from 
said selector if TH2 < r < TH1 ; 

10. The apparatus set forth in claim 8 characterized in that said said apparatus responsive to data representing image 
information of image areas adjacent said block of bad pixel data, for determining a measure, r, of correlation of 
image features, comprises: 

apparatus for determining at least the relative correlation of image gradients in blocks of pixel data adjacent 
said block of bad pixel data. 

1 1 . The apparatus set forth in claim 1 0 characterized in that said said apparatus responsive to data representing image 
information of image areas adjacent said block of bad pixel data, for determining a measure, r, of correlation of 
image features, further comprises: 

apparatus for determining the relative correlation of image motion in blocks of pixel data adjacent said block 
of bad pixel data; and 

apparatus for providing the weighted sum of values representing said relative correlation of image motion 
and said relative correlation of image gradients to generate said correlation measure, r. 

1 2. The apparatus set forth in claim 1 1 characterized in that said said apparatus responsive to data representing image 
information of image areas adjacent said block of bad pixel data, for determining a measure, r, of correlation of 
image features, further comprises: 

means for determining a temporal correlation measure between a block adjacent said block of bad pixel data 
and a block corresponding to a predicted block of a block adjacent said block of bad pixel data; and 

apparatus for providing the weighted sum of values representing said relative correlation of image motion, 
said relative correlation of image gradients and said temporal correlation measure, to generate said correlation 
measure, r. 

1 3. The apparatus set forth in claim 1 1 characterized in that said said apparatus for determining the relative correlation 
of image motion includes apparatus for calculating the relative correlation of motion vectors associated with blocks 
of compressed video data adjacent said block of bad pixel data. 

14. A method for generating a block of pixel data for a missing or corrupt block of pixel data (hereafter bad block), char- 
acterized by; 

providing blocks of pixel data and an error signal indicating which of said blocks are bad blocks; 
determining a dominant image gradient in blocks of pixels adjacent a bad block; 

spatially interpolating, in accordance with said dominant image gradient, a substitute block of data for said 
bad block in response to said error signal; 

temporally predicting a substitute block of data for said bad block in response to said error signal; 
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converting the spatially interpolated substitute block to a block of transform coefficients, TC1, representing 
frequency spectra; 

converting the temporally predicted substitute block to a block of transform coefficients, TC2, representing 
frequency spectra; 

forming a further block of coefficients by selecting respective coefficients from said blocks of transform coef- 
ficients TC1 and TC2 according to a predetermined criteria; 

inverse transforming said further block of coefficients to produce a substitute block of pixel data. 

15. The method set forth in Claim 14, characterized by determining a dominant image gradient in blocks of pixels adja- 
cent a bad block, and spatially interpolating in accordance with said dominant image gradient; wherein said step of 
selecting comprises: 

selecting from said block of transform coefficients TCI, only coefficients representing image information 
along the dominant gradient; and 

selecting from the block of transform coefficients TC2, the remaining coefficients. 

16. The method set forth in claim 14 characterized in that said the step of selecting comprises: 

providing an error value "e"; 

determining differences between corresponding coefficients in the blocks of transform coefficients TCI and 
TC2; and 

if the respective difference is less than the error value "e", selecting the coefficient from the block of trans- 
form coefficients TC1 , else selecting the respective coefficient from the block of transform coefficients TC2. 

1 7. The method set forth in claim 16 characterized in that said the step of providing the value "e", comprises: 

providing a block of data adjacent the bad block, and a temporally predicted block corresponding to said 
block of data adjacent said bad block; 

determining differences between corresponding data in said temporally predicted block and the correspond- 
ing block; and 

statistically analyzing the differences and generating the value "e" from a statistic of the differences. 
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